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ABSTRACT 

We present measurements of the effective temperatures and surface gravities for 
a sample of hot DA white dwarfs, using the Lyman line data available from the HUT, 
ORFEUS and FUSE far-UV space missions. Comparing the results with those from the 
standard Balmer line technique, we find that there is a general good overall agreement 
between the two methods. However, significant differences are found for a number of 
stars, but not always of a consistent nature in that sometimes the Balmer temperature 
exceeds that derived from the Lyman lines and in other instances is lower. We conclude 
that, with the latest model atmosphere calculations, these discrepancies probably do 
not arise from an inadequate theoretical treatment of the Lyman lines but rather from 
systematic effects in the observation and data reduction processes, which dominate 
the statistical errors in these spectra. If these systematic data reduction effects can 
be adequately controlled, the Lyman line temperature and gravity measurements are 
consistent with those obtained from the Balmer lines when allowance is made for 
reasonable observational uncertainties. 

Key words: stars:atmospheres - stars:white dwarfs - ultraviolet:stars. 



* 



1 INTRODUCTION 

Two of the most fundamental measurements required to un- 
derstand the nature of any star are the determination of 
its surface gravity and effective temperature. In the case of 
the white dwarf stars, where nuclear burning has long since 
ceased, temperature is an indication of their cooling age and, 
therefore, enables us to map out the evolutionary sequence. 
Early studies of the white dwarf population relied on pho- 
tometric measurements (e.g. Koester et al. 1979). However, 
a major break through in the reliability and accuracy of 
these measurements was achieved with the development of 
a spectroscopic technique for the H-rich DA white dwarfs, 
pioneered by Holberg et al. (1986), where the IUE observed 
Hi Lyman a lines of a number of DA white dwarfs were 
compared with synthetic model stellar atmospheres. Similar 
techniques can be applied to the H i Balmer lines (Holberg 
et al. 1985). Indeed, when multiple Balmer lines (usually 3 
or 4) are systematically analysed a unique well-defined de- 
termination of both the temperature and surface gravity can 
be obtained. Bergeron et al. (1992, BSL) were the first to 



apply this technique to a large sample of DA white dwarfs. 
Combining the measurements of T B s and log g with the evo- 
lutionary calculations of Wood (1992), BSL obtained esti- 
mates of the mass for each star, confirming the narrowness 
of the white dwarf mass distribution and yielding an accu- 
rate measurement of its peak. Subsequent studies of large 
samples of white dwarfs from both optical and EUV sur- 
veys have firmly established this technique as the primary 
way of determining the DA white dwarf temperature scale 
and placing these objects in their evolutionary context (see 
Vennes et al. 1997; Marsh et al. 1997; Finley et al. 1997a). 
The Balmer line technique can only be applied if these 
features are visible in the stellar spectrum. However, there 
are several situations where this is not the case or where the 
measurement is compromised in some way. For example, if 
a white dwarf resides in a binary system with a much more 
luminous main sequence or evolved companion, the optical 
spectrum will be dominated by the latter star. A well-known 
illustration of this is the DA+K star binary V471 Tauri, 
which has been studied extensively (e.g. Barstow et al. 1997 
and references therein) and where the Balmer lines are not 
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Figure 1. Section of the HUT spectrum of G191— B2B, covering 
the wavelength range from 900A to 1200A. 



detectable in the glare of the K2 dwarf. One major result 
of the EUV sky surveys has been the discovery of many 
more similar systems with companion spectral types ranging 
from A to K (e.g. Barstow et al. 1994a; Vennes et al. 1998; 
Burleigh et al. 1997). 

While the signature of the white dwarf cannot be vi- 
sually separated from its companion in these binaries, the 
UV flux from the degenerate star dominates that from the 
late-type star, provided the primary has a spectral type later 
than ~ A2. Consequently, most of the white dwarfs in such 
systems have been identified from IUE spectra. In principle, 
the white dwarf temperature and gravity can be determined 
from the Lyman a line profile in each case. However, a sin- 
gle line is unable to give an unambiguous measurement of 
both T e ff and log;; simultaneously. Additional information, 
such as the stellar distance, can constraint the possible val- 
ues further, but the distance measurements may not always 
be sufficiently accurate or may be based on uncertain knowl- 
edge of the primary spectral type. The importance of access 
to the Lyman series in an unresolved binary is illustrated 
by determination of T c g and logg for several important 
white dwarfs including V471 Tauri (Barstow et al. 1997) and 
HZ43 (Dupuis et al. 1998), based on spectra obtained by the 
ORFEUS spectrometer (see below). More recently, Burleigh 
et al. (2001) have obtained T c s and log g for the white dwarf 
companion to the A star j3 Crateris, using FUSE. 

In principle, the full Lyman line series could be used in 
the same manner as the Balmer lines to determine T c g and 
log g. While the IUE and HST bands do not extend to short 
enough wavelengths to encompass more than Lyman a, the 
short duration HUT and ORFEUS missions have provided 
several observations of white dwarfs down to the Lyman 
limit with which to examine this idea. The HUT telescope 
was carried into space twice aboard the Space Shuttle, as 
part of the Astro 1 and Astro 2 missions. A total of eight 
isolated hot DA white dwarfs were observed at a resolu- 
tion of ~ 4A (fwhm). Launched from the Space Shuttle, the 
ORFEUS missions utilised a free-flying far-UV spectrome- 
ter operating at a factor 10 higher resolution (~ 0.3A) than 
HUT. A total of four DAs were observed during two separate 
flights. 

Several authors have used these far-UV spectra to de- 
termine T e ff and log g. However, as the total number of stars 
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Figure 2. ORFEUS spectrum of G191-B2B. 

studied on any one mission has been small, it has been dif- 
ficult to establish whether or not the results of Lyman line 
analyses are generally in agreement with Balmer line ob- 
servations of the isolated DA stars. A test of this nature is 
crucial if we are to routinely use Lyman line data for this 
purpose. A preliminary investigation of this issue was carried 
out by Finley et al. (1997b), using the Astro 2 HUT data. 
Their major conclusion was that the Lyman lines were sig- 
nificantly weaker than the standard Stark broadening theory 
of Vidal et al. (1973). Conflicting results have been obtained 
with the ORFEUS data. For example, Barstow et al. (1998), 
utilizing state-of-the-art non-LTE models, found that the 
Lyman and Balmer T c g determinations were in good agree- 
ment for the heavy element rich star G191— B2B but not for 
the similar object REJ0457-281. Dupuis et al. (1998) also 
achieved good agreement between the results of their Lyman 
line analysis of the pure H atmosphere white dwarf HZ43 
and other published values of T e ft from a variety of pub- 
lished sources including HUT, EUVE besides ground-based 
Balmer line observations. However, it must be remembered 
that these measurements reported in the literature were de- 
rived using several different model codes, spanning several 
generations of these programmes. Hence, such comparisons 
lack the uniformity and self-consistency required to evalu- 
ate the efficacy of the Lyman lines as reliable T c g and log g 
indicators. The launch of the Far Ultraviolet Spectroscopic 
Explorer (FUSE) on 1999 June 24 has provided long dura- 
tion access to the Lyman series region of the electromagnetic 
spectrum for the first time since Copernicus, in the early 
1970s. However, Copernicus was not sufficiently sensitive to 
observe any white dwarfs, except Sirius B. With a spectral 
resolution superior to both HUT and ORFEUS, and compa- 
rable effective area, FUSE promises to produce Lyman series 
data for many white dwarfs. This paper presents a critical 
analysis of the use of the H Lyman series absorption lines 
to determine T e ff and log g. We reanalyze the archival HUT 
and ORFEUS spectra in conjunction with new, improved 
non-LTE synthetic spectra, comparing the results with our 
archive of Balmer line data, using the same models. We also 
consider the possible systematic observation and data re- 
duction effects that might affect either or both Balmer and 
Lyman procedures. 
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Table 1. Far-UV Lyman series spectra used in this study, ob- 
tained from HUT, ORFEUS and FUSE missions. 



Star 


Instrument 


Data set name/number 


GD50 


HUT2 


A19001, A19002 


GD394 


HUT2 


A09201, A09202, A13801, A21101 


HZ43 


HUT 


A32901, A32902 




ORFEUS2 


12, 3270659-3362240 


REJ0512 


HUT2 


A12701, N12701 


G191-B2B 


HUT2 


N11401 




ORFEUS1 


4, gl91b2b.l-4 


GD153 


FUSE 


M1010401 


GD71 


HUT2 


N11201 




FUSE 


M1010301, M1010302 


REJ0457 


ORFEUSl 


4, mct0455_l-4 


PG1342+444 


FUSE 


A0340402 


REJ0558 


FUSE 


A0340701 


REJ1738 


FUSE 


A0340301 
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Figure 3 — continued b) After processing, as described in the 
text and resampled with a 0.1 A bin size to improve the s/n. 




900 920 910 960 9B0 1000 1040 1080 

Wavelength (A) 



Figure 3. FUSE spectrum of PG1342+444, a) combined from 
the individual grating spectra, before applying our processing 
pipeline. 



2 OBSERVATIONS 

2.1 Lyman line spectra 

We have obtained far-UV spectra covering the Lyman se- 
ries from the /3 line through to the series limit from three 
telescopes, HUT, ORFEUS and FUSE. Table [j] summarises 
all the observations, which are discussed briefly below. All 
the spec tra were obtained from the Multi- mission archive 
(MAST, http://archive. stsci.edu/mast.htmj ) hosted by the 
Space Telescope Science Institute, or directly from the in- 
vestigators of the respective instruments. 



2.1.1 HUT observations 

The HUT instrument was flown on two space shuttle mis- 
sions, Astro 1 and Astro 2, in 1990 December and 1995 
March respectively. The payload consisted of a 90cm f/2 tele- 
scope with a Rowland circle spectrograph at the prime focus, 
covering the 820-1840A wavelength range, in first order. The 
spectral resolution (A/AA = 500) was dependent upon the 
instrument temperature and pointing stability. The final cal- 
ibration is described by Kruk et al. (1997) The second flight 
incorporated significant improvements in sensitivity and sta- 
bility. The spectral resolution was wavelength dependent in 
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Figure 4. Ground- based optical spectrum of GD71. 



this flight, varying from 2-4.5 A (fwhm). Both temporal and 
wavelength dependent sensitivity changes were experienced 
during the mission, which were monitored and character- 
ized by multiple observations of three white dwarfs. Kruk 
et al. (1999) discuss the final Astro 2 HUT calibration. An 
example HUT spectrum, of G191-B2B, is shown in figure |j 
While eight stars were observed, in RE J 1738 strong H2 ab- 
sorption prevents sensible analysis of the H Lyman lines (the 
increased resolution allows the H2 lines to be removed in the 
FUSE spectrum). In addition, Wolf 1346 is too cool to show 
a full Lyman line series and is not suitable for this study. 
Hence, we only made use of observations of six stars. 



2.1.2 ORFEUS observations 

Higher resolution far-UV spectra were obtained with the 
ORFEUS telescope (Hurwitz & Bowyer 1991) on the free- 
flying Astro-SPAS platform, deployed from the shuttle. OR- 
FEUS was flown on two occasions, a 5 day mission in 1993 
September and a 14 day mission in 1996 November. On the 
first flight, the Berkeley spectrometer located at the focus 
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of the 1-m primary, covered the spectral range 360-1176A 
at a resolution of A/AA = 5000 (Hurwitz & Bowyer 1995; 
Hurwitz et al. 1998). For the second mission, the far-UV 
wavelength coverage was 900-1200A and the spectral resolv- 
ing power 3300 (Hurwitz et al. 1998). 

All data obtained during both ORFEUS missions are 
now in the public domain. The raw spectra suffer from con- 
tamination by a scattered light component, which is a com- 
bination of EUV flux from the second spectral order and 
direct scatter from the grating (Hurwitz, private commu- 
nication). When added to the stellar spectrum, the appar- 
ent level of the continuum flux will be higher than the true 
one. In determinations of T e ff and \ogg, the line profiles, 
whether Balmer or Lyman, are measured from the contin- 
uum level. Therefore, the scattered light component must be 
accounted for in any analysis to avoid obtaining erroneously 
high temperatures from apparently weaker lines. Immedi- 
ately shortward of the Lyman limit, the stellar spectrum 
makes no contribution to the net flux, as a result of absorp- 
tion by interstellar neutral hydrogen. Therefore, the level of 
the scattered light component was estimated from the ob- 
served flux in the 850-900A range. The ORFEUS spectrum 
of G191— B2B is shown in figure H, for comparison with the 
HUT spectrum (figure hi). While four stars were observed, 
one of these, the white dwarf in V471 Tauri is in a close bi- 
nary and no Balmer series data can be obtained due to the 
overwhelming brightness of the K star companion. 

2.1.3 FUSE spectra 

The FUSE mission was placed in low Earth orbit on 1999 
June 24. After several months of in-orbit checkout and cal- 
ibration activities, science observations began during 1999 
December. An overview of the FUSE mission has been given 
by Moos et al. (2000) and the spectrograph is described in 
detail by Green et al. (1994). Further useful information is 
included in the FUSE Observer's Guide (Oegerle et al. 1998) 
which can be found with other technical documentation on 



Table 2. Nominal wavelength ranges (A) for the FUSE detector 
segments 



the FUSE website (ittp://fuse. pha.jhu.edu). 

Only a few scientific papers have, as yet, been published 
incorporating FUSE data. Hence, it is appropriate to give 
a brief description of the spectrometer and its current sta- 
tus in the context of the reduction and analysis of the data 
presented here. The far-UV spectrometer is based on the 
Rowland circle design and comprises four separate optical 
paths (or channels) . Each channel consists of a mirror, a fo- 
cal plane assembly (including the spectrograph apertures), 
a diffraction grating and a section of one of two detectors. 
To maximise the throughput of the instrument, the chan- 
nels must be co-aligned so that light from a single target 
properly illuminates all four channels. Two of the mirrors 
and two of the gratings are coated with SiC to provide sen- 
sitivity at wavelengths below » 1020A, while the other two 
mirror/grating pairs are coated with Al and a LiF overcoat. 
This latter combination yields about twice the reflectivity 
of SiC at wavelengths above 1050A, but has little reflectiv- 
ity below 1020A. The overall wavelength coverage runs from 
905A to 1187A. 

Spectra from the four channels are recorded on two mi- 
crochannel plate detectors, with a SiC spectrum and LiF 
spectrum on each. The individual detectors are divided into 
two functionally independent segments (A and B) , separated 



Channel Segment A 

SiC 1 1091.1-1003.9 

LiF 1 987.1-1082.2 

SiC 2 916.8-1006.4 

LiF 2 1181.7-1085.6 



Segment B 

992.6-905.0 

1094.3-1187.7 

1016.3-1105.0 

1074.6-978.1 



by a small gap. Consequently, there are eight detector seg- 
ment/spectrometer channel combinations to be dealt with 
in reducing the FUSE data. The nominal wavelength ranges 
of these are listed in table bl 

Several problems have been reported during in-orbit op- 
erations of the FUSE satellite which need to be taken into 
account in the data reduction process. Maintaining the co- 
alignment of the individual channels has been difficult, prob- 
ably due to thermal effects. Sometimes, during an observa- 
tion, a target may completely miss the aperture of one or 
more channels, while being well centred in the others. In 
addition, even if the channels are well-aligned at the begin- 
ning of an observation, the target may subsequently drift 
out of any of the apertures. To minimize this problem, most 
observations have been conducted using the largest aper- 
ture available (LWRS, 30 x 30 arcsec), limiting the spec- 
tral resolution for point sources to between 10000 and 20000 
(« 0.05 to 0.1 A) compared to the 24000-30000 expected for 
the 1.25 x 20 arcsec HIRS aperture. 

Most of the FUSE spectra included in this work were ob- 
tained in time-tag mode (where the arrival time of each de- 
tected photon is recorded and the image subsequently recon- 
structed from the positional data also included in the data 
stream) in several separate exposures through the LWRS 
aperture (table hi). However, two observations (GD71 and 
GD153) were obtained from the MAST archive as part of the 
early release of data for calibration purposes. All these data 
were processed with the CALFUSE pipeline version 1.6.6. 
Subsequently, a new wavelength calibration has been made 
available and which is applied to all data processed recently 
(pipeline version 1.7). We have replaced the original wave- 
length calibration files with the revised versions according 
to the procedure described on the FUSE data analysis web 
pages (fuse.pha.jhu.edu/analysis/wavelength_062200.html), 
before combining the individual spectra. 

Initially, for each star, we considered the separate 
exposures for a single channel/detector segment. Since, 
the signal-to-noise of these is relatively poor and the 
wavelength binning (~ 0.006A) over samples the true 
resolution by a factor of 2-3, all the spectra were re- 
binned to 0.02A pixels for examination. Fortuitously, sev- 
eral prominent interstellar absorption lines are detected 
in each spectrum showing that, for a given channel and 
detector side, the wavelength scales of each exposure are 
well aligned. Consequently, it is a straightforward process 
to co-add the individual exposures, to produce a single 
spectrum. We used the iraf script FUSECOMBINE (see 
fuse.pha.jhu.edu/analysis/IRAF_scripts.html) to co-add the 
multiple exposures, which weights the individual spectra by 
their exposure time. This whole procedure was repeated for 
all eight channel/detector side combinations for each star. 

Since the complete FUSE wavelength range is covered 
in 3 ~ IOOA bands, a number of the individual detector 
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segments overlap almost completely in wavelength. For ex- 
ample, the SiC 1A, SiC 2B, LiF 2A and LiF 2B all cover the 
range « 1000-llOOA, while the SiC 2A/2B (w 900-1000A) 
and LiF 1A/1B (« 1100 - 1200A) match in wavelength. 
Consequently, to achieve the optimum signal-to-noise for 
analysis of the data, it is desirable to combine all these 
individual segment spectra. We have written a small For- 
tran programme to do this, which is able to take into ac- 
count the differing wavelength ranges and spectral bin sizes 
of each. First, all the spectra are re-sampled onto a com- 
mon wavelength scale and then re-binned into 0.02A steps 
to reduce the level of over sampling and avoid any fringing 
effects that might arise from the first part of the proce- 
dure. The re-sampled/re-binned spectra are then co-added, 
weighting individual data points by the statistical variance, 
averaged over a 10A interval, to take into account the differ- 
ences in effective area of each segment and any differences 
in exposure time that may have arisen from rejection of bad 
data segments. We find, through visual inspection of the 
observed spectra, that the statistical noise tends to increase 
towards the edges of the wavelength range. In cases where 
the signal-to-noise is particularly poor in these regions, we 
have trimmed the spectra to remove these data points prior 
to co-addition. 

As an example, the FUSE spectrum of PG1342+444 is 
shown in figure Ha, combined from the individual grating 
spectra but without any of the processing described above. 
The strong high flux spikes are geocoronal emission lines, 
while poor s/n regions of the spectrum (usually covered by 
only one grating) are revealed by increased scatter in the 
data points (e.g. near 1080 A). A few of the stronger heavy 
element and interstellar absorption lines are also visible. The 
broad dip in the 1160- 1180 A region of the spectrum is an 
artefact known as the "worm", which is a local 10-20% re- 
duction in the effective area and produces a corresponding 
loss of flux. Figure yb shows the same spectrum processed 
according to our prescription and resampled with a bin size 
of 0.1 A to improve the s/n for our analysis. 

2.2 Balmer line spectra 

The majority of the optical data we use here was ob- 
tained as part of a spectroscopic follow-up programme fol- 
lowing the ROSAT X-ray and EUV sky survey. Observa- 
tions were undertaken in both Northern and Southern hemi- 
spheres. Southern hemisphere data were obtained with the 
1.9-m Radcliffe reflector of the South African Astronomi- 
cal Observatory (SAAO), while stars in the Northern Hemi- 
sphere were observed with the Steward Observatory 2.3-m 
telescope on Kitt Peak. Full details of these observations 
have already been published by Marsh et ai. (1997). The 
main difference between the Southern and Northern Hemi- 
sphere data is their spectral resolution, « 3A (FWHM) and 
w 8A (FWHM) respectively. Some of our original Northern 
Hemisphere spectra did not cover the complete Balmer line 
series, excluding H/3, due to the limited size of the CCD chip 
available at the time. Hence, in those cases, we have replaced 
the archival data with more recent observations (GD394 
and GD71) made using the same instrument but with a 
larger chip (see e.g. figure ||). For some stars (GD50, HZ43, 
G191-B2B, and GD153), while archival ("old") data are of 



«l 



u 



o 

X 

CN 

m 



E" o 

CD ' — 
X 

x m 



V 




\r*~ 



4000 



4500 
Wavelength (A) 



Figure 5. Fit to the ground-based optical H/3 to He lines of 
G191-B2B (T efi = 51510 ± 880K, logg = 7.53 ± 0.09). 
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Figure 5 - continued b) Fit to the HUT L/3 to Le lines of 
G191-B2B (T cH = 52930 ± 360K, logs = 7.16 ± 0.20). Although 
the line cores are not included in the analysis, they are shown 
here for illustrative purposes. 



good quality and cover all the Balmer line series the "new" 
observations have longer exposures and, consequently, better 
s/n. We analyse all these data, as they provide a useful test 
of consistency between repeated observations of the same 
star, using basically the same instrument configuration but 
on different nights. 



3 MODEL ATMOSPHERE CALCULATIONS 

One potential flaw in comparing the results from published 
Lyman and Balmer line analyses is that different authors 
may have utilized different stellar atmosphere codes, even 
spanning successive generations of these codes. To ensure 
that our work is at least internally consistent we re-examine 
all the Balmer line data with the latest version of the stellar 
atmosphere programme tlusty and its associated spectral 
synthesis package synspec . 
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Figure 5 - continued c) Fit to the ORFEUS L/3 to he lines of 
G191-B2B (T cff = 53180 ± 520K, logg = 7 A3 ± 0.04). Although 
the line cores are not included in the analysis, they arc shown 
here for illustrative purposes. 



We have calculated completely new grids of model stel- 
lar atmospheres using the non-LTE code tlusty (Hubeny & 
Lanz 1995). These are based on work reported by Lanz et al. 
(1996) and Barstow et al. (1998, 1999). Two separate sets of 
calculations were performed: pure H models, for those stars 
without significant abundances of heavy elements, and mod- 
els with a homogeneous mixture of heavier elements, includ- 
ing C, N, O, Si, Fe and Ni, for the others. The upper limit to 
the temperature of the pure H models was 70000K, to span 
the range within which DAs with pure H atmospheres are 
found, while the heavy element calculations were extended 
to 120000K, to deal with hotter DA stars, although we note 
that in this particular study the hottest star we consider is 
REJ1738, with T off « 70000K. To take account of the higher 
element ionization stages that are likely to be encountered 
in the latter objects, we have added new ions of OVI, Fe- 
VII/VIII and NiVII/VIII to the model atoms as well as ex- 
tending the data for important ions such as CIV to include 
more energy levels. As before, all the calculations were per- 
formed in non-LTE with full line-blanketing, including Stark 
broadening of all the CNO lines. 

The stars included in this study are divided into two 
distinct groups. Those with pure H envelopes, for which 
we used the pure H model calculations, and those with sig- 
nificant quantities of heavy elements. For the latter group, 
we fixed the abundances of the heavy elements at the val- 
ues determined from our earlier homogeneous analysis of 
G191-B2B (He/H=1.0xlO- 5 ,C/H=4.0xl(r 7 ,N/H=1.6x 
1(T 7 , 0/H=9.6 x 1(T 7 , Si/H=3.0 x 1(T 7 , Fe/H=1.0 x 1(T 5 , 
Ni/H=5.0 x 10~ 7 ), but taking into account that the CIV 
lines near 1550A have subsequently been resolved into mul- 
tiple components by STIS (Bruhweiler et al. 2000; Bannis- 
ter et al. 2001). While not all stars have exactly the same 
heavy element abundances, our recent work has shown that 
the differences are not very large (see Barstow et al. 2001) 
and, at this level, will not have a significant effect on the 
Balmer/Lyman line measurements (see e.g Barstow et al. 
1998). 



In the spectrum synthesis code synspec (Hubeny et al. 
(1994), we have replaced the hydrogen Stark line broaden- 
ing tables of Schoning & Butler (private communication) 
by the more extended tables of Lemke (1997). The latter 
allow a more accurate interpolation of the electron density 
for high density environments, such as the atmospheres of 
white dwarfs. The spectra produced by the tlusty /synspec 
codes were recently extensively tested against the results of 
Koester's codes (Hubeny & Koester in preparation). The 
differences in the predicted spectra for T e g =60000K and 
log g — 8 were found to be below 0.5% in the whole UV and 
optical range. Furthermore, we have found that the inaccu- 
racy in the interpolations of the Schoning and Butler tables, 
together with some fine details of our treatment of the level 
dissolution, were the primary reason for the disagreement 
between the spectroscopically deduced T e g using tlusty and 
Koester models obtained by Barstow et al. (1998). These 
changes largely resolve the differences between codes noted 
by Bohlin (2000). 



4 DETERMINATION OF TEMPERATURE 
AND GRAVITY 

The technique for determining T e s and logg, by comparing 
the line profiles with the predictions of synthetic spectra 
is well established (see Holberg et al. 1986; Bergeron et al. 
1992 and many subsequent authors). We have described our 
own Balmer line analysis technique in several earlier papers 
(e.g. Barstow et al. 1994b), but as the results presented in 
this paper rely heavily on it we repeat the details here. The 
same technique can also be applied to analysis of the Lyman 
lines, as has already been demonstrated (e.g. Barstow et al. 
1997; Barstow et al. 1998). However, we have modified our 
earlier approach, as outlined below. 

Both sets of T e g and logg measurements were per- 
formed using the program xspec (Shafer et al. 1991), which 
adopts a \ 2 minimization technique to determine the model 
spectrum which gives the best agreement with the data. For 
the Balmer lines, the four strongest (/3, 7, 8, e) are simultane- 
ously included in the fit and an independent normalisation 
constant applied to each, ensuring that the result was inde- 
pendent of the local slope of the continuum and reducing the 
effect of any systematic errors in the flux calibration of the 
spectrum. Individual lines maintain their local continuum 
slope but are decoupled from the overall energy distribution 
of the entire Balmer line region. 

In the case of the Lyman lines, the analysis needs to be 
handled in a slightly different way. In general, at given values 
of T e ff and logg, the Lyman lines are stronger, and, apart 
from the a and /3 lines, overlap substantially. To deal with 
the Lyman data for each instrument, we separated out the 
ft (1000- 1050A) line and incorporated the remaining lines 
(7 through e inclusive) into a single file. To take account of 
any low frequency systematic effects in the flux calibration 
we applied individual normalization constants to each of the 
two sections of data. 

Since the Lyman data are, by definition, obtained from 
space-based platforms, there are no systematic errors arising 
from an atmospheric extinction correction. Furthermore, the 
flux calibration is usually obtained from a detailed off-line 
calibration of the instrument, applied as part of a standard 
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Figure 6. Scatter plot of the values of T e g measured using the Balmer and Lyman series lines. Diamonds are HUT data, the squares 
are the ORFEUS data and the triangles the FUSE data. The error bars displayed correspond just to the statistical la uncertainties. The 
upper solid line corresponds to equal Balmer and Lyman temperatures, while the lower line is the least squares straight line fit to the 
data, corresponding to y=0.964x. 



pipeline, rather than direct comparison of the observed spec- 
trum with that of a selected standard star. However, the 
instrument calibration still refers to observations of well- 
studied stars, such as white dwarfs, some systematic uncer- 
tainties will still apply, but will be different to those arising 
from the ground-based techniques Although, in general, no 
extinction correction needs to be applied to the Lyman se- 
ries data, there are two important effects that need to be 
accounted for in any analysis. First, the cores of the Lyman 
absorption lines can be contaminated by geocoronal emission 
components. Second, interstellar absorption can artificially 
deepen the core of the stellar Lyman absorption lines. These 
two effects compete with each other and may occasionally 
conspire to cancel each other out, but usually they must be 
removed from the data in an appropriate way. 

The observed strengths of individual geocoronal lines 
depend on a number of factors, including the strength of 
the Solar line being scattered, the density of scattering 
atoms along the line of site (here Hi atoms) and the opti- 
cal depth of the scattering environment. In general however, 



the strongest line is always Lyman a followed by Lyman /3. 
In most of the data considered here, the Lyman a line is 
not observed and is only available in the HUT data. Hence, 
it is not used for consistency of analysis and the Lyman j3 
line is by far the strongest component we need to consider 
(e.g. Figure H) . From a signal-to-noise point of view the rel- 
ative strength of the emission when compared to the stellar 
continuum is also affected by the viewing geometry, the in- 
strument design and the brightness of the target star. The 
geocoronal radiation has a natural spectral width caused by 
Doppler effects, which is then further broadened in the in- 
strumental aperture by the diffuse nature of the source and 
the spectral resolution of the instrument. As a result of pos- 
sible differences in the relative velocity of the stellar and 
geocoronal sources, the emission lines are not necessarily 
aligned with the stellar absorption cores. 

To make sensible estimates of the uncertainty in the fit- 
ted parameters, the value of the reduced \ 2 should be less 
than w 2. This is the case for the fits to all the spectra in- 
cluded in this analysis, which are, therefore, formally good 
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Figure 7. Scatter plot of the values of T e g measured using the "old" (x-axis) and "new" (y-axis) Balmer series lines. The error bars 
displayed correspond to statistical lc uncertainties alone. The light solid line corresponds to equal Balmer temperatures. 



fits. Errors on T c g and log g can then determined by allow- 
ing the model parameters to vary until the 5\ 2 reached the 
value corresponding to the la level for 2 degrees of freedom 
(2.3). It should be noted that these estimates only include 
statistical uncertainties and do not take into account any 
possible systematic effects related to the data acquisition 
and reduction processes. 

Figure ]a shows examples of the fits to the Balmer 
(5a) and Lyman (5b, 5c) lines of G191— B2B respectively. 
G191— B2B is one of only two stars (the other is HZ43), 
where spectra exist for ground-based, HUT and ORFEUS, 
providing a useful cross-check between all these instruments. 
This is particularly important as, prior to the launch of 
FUSE, the majority of the Lyman line DA data had been 
obtained by HUT. 

Table H summarises the effective temperatures and 
gravities obtained from the spectral fits for all data sets for 
each star, including the 1 a statistical uncertainties on these 
values. Where more than one ground-based balmer line mea- 
surement is reported, the upper value is obtained with the 



"old" archival data and the lower is from the "new" , higher 
s/n observation. For GD394 and GD71, the reported mea- 
surements are only from the "new" spectra as the archival 
obervations only span 2-3 lines of the Balmer series and are 
not used here. 

The T ff and logg measurements obtained from the 
various observations are a somewhat heterogeneous group. 
While we have a complete set of ground-based optical spec- 
tra for all the stars, and in some cases multiple ground-based 
observations of an individual star, the far-UV observations 
are incomplete, with different instruments covering a dif- 
ferent group. At this stage, for no star is there a complete 
set of observations with data from all space-based instru- 
ments. In a number of cases, there are multiple observations 
of a star made with a single telescope. Rather than combine 
these spectra into a single average, where the signal-to-noise 
is adequate, we have chosen to analyse these individually to 
see if any systematic differences may occur between separate 
observations of the same star with the same instrument. We 
discuss the comparison of individual groups of data below. 
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Table 3. Best-fit effective temperature and surface gravity for each observation of each DA white dwarf in the sample. Multiple values 
given each column refer to independent observations using the same instrument. 



STAR GROUND HUT 

BASED 

T cS (K) (err) log g (err) T cff (K) (err) log g (err) 

GD50 39060(950) 9.30(0.12) 37300(600) 8.81(0.19) 

42040(570) 9.15(0.05) 38200(600) 8.82(0.18) 

GD394 39290(360) 7.89(0.05) 34790(300) 7.96(0.10) 

33100(450) 7.63(0.15) 

34200(500) 8.01(0.10) 

34150(400) 8.00(0.10) 

HZ43 50370(780) 7.85(0.07) 46100(1300) 7.65(0.09) 

48500(400) 8.05(0.02) 47800(1250) 7.84(0.12) 

45200(1000) 7.68(0.13) 

RE0512 31670(140) 7.20(0.04) 31100(250) 7.21(0.10) 

31100(350) 7.19(0.14) 

G191B2B 51510(880) 7.53(0.09) 52930(3600) 7.16(0.20) 

54560(200) 7.60(0.02) 

GD153 39290(340) 7.77(0.05) 

37120(140) 8.02(0.02) 

GD71 32780(65) 7.83(0.02) 31970(150) 7.72(0.07) 



REJ0457 


50960(1070) 


7.93(0.08) 


PG1342 


66750(2450) 


7.93(0.11) 


REJ0558 


59510(2210) 


7.70(0.14) 


REJ1738 


66760(1230) 


7.77(0.10) 



OREErjS(o)/Ft/SE(f) 
T cB (K) (err ) 



(o)50080(290) 



logg (err) 



.14(0.05) 



(o)53180(520) 


7.43(0.04) 


(f)37150(100) 


8.1(0.03) 


(f)32820(1280) 
(f)32220(310) 


8.85(0.28) 
8.49(0.07) 


(o)59300(1150) 


7.57(0.08) 


(f)54286(1000) 


7.82(0.07) 


(f)61140(660) 


7.61(0.05) 


(f)70060(800) 


8.00(0.01) 



5 DISCUSSION 

5.1 Effective temperatures 

In figure o, we show a scatter plot of the values of T e g 
measured by Balmer line (x-axis) and Lyman line (y-axis) 
techniques. All combinations of ground and far-UV obser- 
vations are shown, to illustrate the full range of possible 
systematic differences. Also shown is the straight line corre- 
sponding to equal Balmer and Lyman temperature together 
with the best fit straight line corresponding to y=0.964x. In 
general, there is good agreement between the two measure- 
ment techniques, but there also exist some clear anomalies. 
For example, the Lyman line values of T e s for GD394 are 
consistent with each other, but systematically lower than 
the Balmer line value (T cff balmer = 39290) by w 5000K. 
Similarly, the value of T B g measured from the Balmer lines 
of PG1342+444 is 12000K higher than the corresponding 
Lyman line measurement. However, such differences are not 
all in the same direction. The Lyman line temperature of 
REJ0457 is some 8000K higher than the Balmer tempera- 
ture. In general, the Lyman T B a measurements are about 4% 
below the Balmer line values. At 50000K, this corresponds 
to a difference of 2000K, which is of a similar magnitude to 
many observational errors. 

The information included in figure r\ does not discrim- 
inate between the source of optical observations used for 



comparions with the far-UV measurements. Although the 
number of repeated optical observations is small, it is still 
relevant to look for systematic differences to assess the mag- 
nitude of systematic uncertainties in the determination of 
T ff . Figure M compares the Balmer line temperatures mea- 
sured using "old" spectra from our archive (x-axis) with 
those of improved s/n from the more recent ( "new" ) obser- 
vations. It is clear that, compared to the size of the formal er- 
rors, the two sets of data, give different results. Observations 
for the same stars analysed with the same models disagree 
for all four stars. The differences range from 5-10%, but are 
not all in the same direction. In fact, the best fit straight 
line through the data points is y=1.01x, only 1% differ- 
ent from the ideal y=x. This result would suggest that the 
ground-based temperature measurements can have system- 
atic errors at the 5-10% level, which dominate the statistical 
errors in spectra of the signal-to-noise used here. Since T e s 
and loggr are correlated, similar systematic effects exist for 
the logg estimates. This is not a new result and Bergeron 
et al. (1992) use multiple observations of the same star to 
quantify the systematic errors. The problem is illustrated in 
figure 0, which shows a summary of spectroscopic observa- 
tions of GD153. This suggests that, for ground-based T c g 
and log g measurements, the statistical uncertainties under- 
estimate the true error by a factor of 2 to 3. Nevertheless, 
the statistical errors do provide a measure of the spectral 
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Figure 8. Published independent ground-based Balmer line measurements of T e g and log g for the white dwarf GD153. 



s/n, the goodness of fit and the sensitivity of the results 
to the line strength and should be reported. However, it is 
not appropriate to take them too literally, particularly when 
propagating them through calculations such as mass and ra- 
dius determinations. 

Figure |6| shows values of T e g derived from the various 
Lyman line analyses and each far-UV telescope is denoted 
by a distinctive symbol (diamonds - HUT, squares - OR- 
FEUS, triangles - FUSE). It is possible to identify particular 
trends associated with each instrument. Most of the FUSE 
T e g values are in reasonable agreement with the Balmer line 
measurements. The only dramatic departure is PG1342+444 
which, as noted earlier, has a 12000K lower Lyman line tem- 
perature. There are too few ORFEUS data points, over too 
narrow a temperature range, to say much about the trends 
in this instrument. Both HZ43 and G191— B2B are in agree- 
ment with the Balmer measurements while the Lyman tem- 
perature of REJ0457-281 is 8000K too high. However, we 
have noted in our earlier work (Barstow et al. 1998), that 
the separate Balmer and Lyman line fits for this star are 
not, from the point of view of the statistical errors and the 
goodness of fit (as determined by the F-test described in 
Barstow et al. 1998), significantly different from the average 
fit to both data sets simultaneously. If we a try a combined 
fit for PG1342, we do not get the same result: the Lyman, 
Balmer and combined fits are all significantly different us- 
ing the same statistical criteria. As discussed earlier in this 
section, the statistical errors are not a very good indication 



of the true uncertainty when systematic effects are consid- 
ered. If we adopt a factor 3 scaling of the statistical errors 
as representative of the overall analysis process, the level of 
disagreement between Lyman and Balmer analyses is within 

2(7. 

Considering the whole sample of observations, the HUT 
data seem to yield Lyman line temperatures that are sys- 
tematically lower than the Balmer line values. This is most 
marked in GD394, but is also seen in HZ43 and GD50. 
Dupuis et al. (2000) have already noted this in their de- 
tailed study of GD394. Although the discrepancies observed 
in the other two stars are smaller, it may not be possible to 
rule out a unique systematic problem with the GD394 data. 
Hence, the observed temperature effect might not be real. 
The lower resolution of the HUT data, when compared to 
the other instruments pose a particular problem, as outlined 
earlier, in terms of taking account of geocoronal contamina- 
tion and interstellar absorption. In addition, removing the 
entire line cores to eliminate these effects does appear to 
have a systematic influence on the outcome of the Lyman 
line analysis. 



5.2 Surface gravity measurements 

Most of the emphasis of this analysis and discussion has 
been on the comparison of temperature measurements. How- 
ever, in determining the evolutionary status and important 
parameters such as mass and radius, the measurement of 
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Figure 9. Scatter plot of the values of logg measured using the ground-based Balmer (x-axis) and far- UV Lyman (y-axis) series lines. 
The diamonds are HUT data, the squares are the ORFEUS data and the diamonds the FUSE data. The error bars displayed correspond 
to statistical la uncertainties. The light solid line corresponds to equal Balmer and Lyman line gravities. 



surface gravity is of equal importance. Figure tn shows the 
measured surface gravities corresponding to the tempera- 
ture data in figure pi The fractional statistical errors on 
logg are larger than for the temperatures and most error 
bars overlap the line of equally Balmer and Lyman line grav- 
ity. However, there are some marked departures from this. 
The surface gravity of GD71 obtained from the Lyman line 
analysis is approximately l.Odex greater than the Balmer 
line result, which does not seem reasonable for such a well- 
studied star. On the other hand the surface gravity measured 
for G191— B2B is much lower (by 0.37dex) than the Balmer 
line result. Like the differences in T c a values, these large dif- 
ferences in log g probably arise from systematic effects in the 
data. However, there are no apparent trends for any partic- 
ular instrument and, as in the temperature measurements, 
no overall systematic departures from the line of equal log g. 



5.3 Systematic effects in the analyis of T c s and 

log g 

One of the main propositions of this paper is that measure- 
ments of Tefr and log g using either Balmer or Lyman series 



lines are prone to systematic errors that are not usually well 
quantified. For the Lyman lines, there is the added com- 
plication of geocoronal emission and interstellar absorption 
modifying the line cores, besides any instrumental effects 
such as scattered light/background contamination. Conse- 
quently, we have simulated some of these effects for selected 
HUT and ORFEUS data sets to examine the possible contri- 
bution to the measurement uncertainties. Fluxes were scaled 
to analyse the effect of errors in effective area calibrations 
and constant fluxes were added/subtracted to study the sen- 
sitivity to accurate background subtraction/scattered light 
corrections. For the HUT simulation, we also carried out an 
analysis including the unresolved line cores. In each case, 
measurements were carried out using the spectral analysis 
techniques described in section for the real data. Table 
summarises the results, which are also shown in figure [id. 

It is clear that there are significant differences between 
the measurements of T e g and log g for each of the simulated 
datasets. The systematic effects that have been introduced 
are typically w 5% of the mean flux. These translate into 
observational errors ~ 2 — 5% for Tgg and ~ 0.5 — 1.0% 
for log p. While the precision of the HUT measurements is 
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Figure 10. Scatter plot of the values of T c s and log g measured for the simulated HUT (diamonds) and ORFEUS (squares) data sets. 
Each data point is labelled with a letter as indicated in table Ul 



Table 4. Values of T e g and log g measured for the simulated HUT and ORFEUS data sets. Letter labels correspond to the data points 
in figure lid 



T eff (K) (err) 

57000(2400) 

58300(2500) 

60220(2300) 

60700(2050) 

61100(950) 

T cff (K) (err) 

51300(300) 

52800(300) 

52800(300) 

54400(250) 



HUT simulation 

logg (err) 

7.41(0.20) 

7.51(0.15) 

7.59(0.14) 

7.58(0.12) 

7.61(0.08) 

ORFEUS simulation 

logg (err) 

7.43(0.03) 

7.47(0.03) 

7.46(0.02) 

7.50(0.02) 



index/comment 

a) basic data 

b) data — constant 5% of mean flux 

c) data scaled by factor 1.05 

d) data + constant 5% of mean flux 

e) basic data, but line cores included 

index/comment 

a) data — constant 5% of mean flux 

b) basic data 

c) data scaled by factor 1.05 

d) data + constant 5% of mean flux 



similar to the level of the systematic errors, the systematic 
effects clearly dominate the statistical uncertainties in the 
ORFEUS analyses. It is interesting to note that the total 
spread of temperature measurements (~ 10000K) for the 
simulated data set is as large as some of the perceived dis- 
crepancies in the real observations. 



6 CONCLUSION 

We have measured the effective temperatures and surface 
gravities for a sample of hot DA white dwarfs using the 
Lyman line data available from the HUT, ORFEUS and 
FUSE far-UV space missions, comparing the results with 
those from the standard Balmer line technique. In general, 
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there is good overall agreement between the two methods. At 
the level of the pure statistical errors, differences are found 
between Lyman and Balmer measurements for a number of 
stars. However, across the sample, the discrepancies are not 
always in the same direction. This is not what would be ex- 
pected if the problems arose from limitations of the stellar 
atmosphere calculations and, in particular, the treatment of 
Lyman and Balmer line broadening. Hence, while the de- 
scription of opacity in the region of the Lyman lines is still 
uncertain, because of uncertainties associated with the oc- 
cupation probability formalism, the description of the line 
broadening of the isolated, lower, Lyman lines seems to be 
in good shape. It is not necessary to make semi-arbitrary ad- 
justments to the input physics, as carried out by Finley et ai. 
(1997b) in their analysis of the HUT data, to force agree- 
ment between the Balmer and Lyman line results. Even so, 
there may well be residual second order effects worthy of 
examination. The best fit line in figure M (excluding prob- 
lematic results noted earlier) hints at a few percent system- 
atic difference between Balmer and Lyman values of T Q s ■ 
However, this is at a similar level to the possible system- 
atic data reduction problems discussed above, which must 
be eliminated before we can consider this further. 

It is well known that optical Balmer line measurements 
can be sensitive to systematic effects in the data reduction 
process, in particular the extinction correction that must be 
applied to take account of atmospheric absorption as well as 
the absolute flux calibration. Although the former problem 
is not present in the far-UV spectra it is apparent that this 
method of temperature and gravity measurement is sensi- 
tive to a number of systematic effects that can compromise 
these results. For example, with the lower spectral resolution 
HUT data, the geocoronal and interstellar contamination of 
the line cores cannot be dealt with very easily, leading to a 
systematic shift in measured values of T c s ■ In a simulation, 
inclusion of the line core regions in the analysis yielded er- 
roneously high values of T e s and logg. In other cases, the 
accuracy of subtraction of instrument background or scat- 
tered light may be important. Simulations of over- or under- 
subtracted background components at the level of 5% of the 
mean flux give shifts of similar magnitude in the measured 
value of T e ff . With high signal-to-noise data available from 
a variety of instrumentation, it is clear that our analyses 
are no longer limited by the statistical errors but by the 
systematic ones arising in the process of data acquisition 
and reduction. Consequently, we conclude that, if these ef- 
fects can be adequately controlled, the Lyman line technique 
gives measurement that are consistent with those obtained 
from the Balmer lines, when allowance is made for realistic 
observational uncertainties. 
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